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Abstract. We report the detection of the elusive carbon- 
chain radical C5N in the dark cloud TMC1 and its tenta- 
tive detection in the circumstellar envelope IRC+10216. 
C5N appears to be two orders of magnitude less abundant 
than the related molecule HC5N and much less abundant 
than expected from current gas phase chemistry models. 
In comparison the HC3N to C3N abundance ratio is of the 
order of 10, in reasonable agreement with model predic- 
tions. 

We have also detected in IRC+10216 two lines arising 
from the C3H radical in its excited 1/4 = 1 state. 
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1. Introduction 

That long carbon chain radicals could be abundant and 
play a large role in interstellar chemistry was first rec- 
ognized with the discoveries of C3N in the circumstel- 
lar envelope IRC+10216 (Guelin & Thaddeus 1977) and 
in the dark cloud TMC1 (Friberg et al. 1980). To date 
seven acetylenic chain radicals, C„H, n = 2 — 8, and five 
cyanopolyyne molecules, HC2 n CN, n = 1 — 5, are identified 
in TMC1 and/or IRC+10216 (Bell et al. 1997, Guelin et 
al. 1997). Surprisingly, no cyanopolyyne radicals heavier 
than C3N were so far detected, despite model predictions 
that at least C5N should be abundant (Herbst et al. 1994). 

The non-detection of C5N was first blamed on a small 
dipole moment and on the lack of spectroscopic data. The 
dipole moment /Xo and the rotation constant B of the 
C-chains radicals depend critically on the nature of their 
electronic ground state. The two lowest states, 2 II and 
2 £, are close in energy (Pauzat et al. 1991). Pauzat and 
co-workers predicted from unrestricted Hartree-Fock cal- 
culations that the 2 n state of C5N lay below the 2 E state 
and that [1q was very small. The line strengths scaling 
with /Xg, it was no wonder that C5N escaped detection. 

Things changed when Botschwina (1996) showed from 
more elaborate coupled cluster calculations that the C5N 
ground state was in fact 2 S and that no was as large as 



3.385 D, even larger than the dipole moment of C3N. It 
then became clear that C5N could be detected at least in 
the laboratory. Kasai et al. (1997) succeeded to synthetize 
this species, to measure its microwave spectrum, and to 
derive its rotational, fine and hyperfine constants, making 
at last possible a sensitive search for C5N in space. In 
this Letter, we report the astronomical detection of this 
radical. 



2. Observations and results 

The long carbon-chain molecules and radicals are nearly 
thermalized in TMC1 and IRC+10216 with rotation tem- 
peratures, T rot , in the ranges 6-10 K and 20-50 K, respec- 
tively. Since, for C 5 N, £fa = 0.067 K, the strongest lines 
in these two sources lie in the A = 1 cm and 3 — 4 mm at- 
mospherical windows. We thus searched for C5N at these 
wavelengths, using the Effelsberg 100-m telescope and the 
Pico Veleta 30-m telescope. 

The Effelsberg observations were made in Octo- 
ber 1997, January and March 1998. The telescope was 
equipped with the K-band maser receiver. The weather 
was mostly clear and the system temperature T sys ~ 100 
K. We used the new AK90 autocorrelator split into two 
20Mhz-wide bands of 4048 channels each. One of the 
bands covered both fine structure components of the 
N = 9 ^ 8 transition of C 5 N (25.250 GHz, Kasai et 
al. 1997), the other the J = 21.5 -> 20.5 transition of 
2 n 3/2 C 8 H (25.227 GHz). The local oscillator frequency 
was switched by ±200 kHz or ±500 kHz and the spec- 
tra folded accordingly. The average spectrum obtained in 
TMC1 at the position of the cyanopolyyne peak (1950.0: 
4 h 38 m 38.6 s , 25°35'45"), smoothed to 15 kHz (0.17 kms" 1 ), 
has a r.m.s. noise of 3.5 mK (units of Tmb)- 

In addition to C5N and CsH, we observed briefly the 
J = 10 -> 9 line of HC 5 N (23.96 GHz) and the N = 2 -> 1 
lines of C 3 N (19.79 GHz). All the data were calibrated 
following the procedure described by Schilke & Walmsley 
(1991) by observing at every frequency the planetary neb- 
ula NGC 7027, whose flux was taken equal to 5.8 Jy. The 
calibration uncertainty, which results mostly from the at- 
mospheric absorption correction and from beam efficiency 
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Fig. 1. Spectrum observed with the 100-m Effels- 
berg telescope toward TMC 1 (cyanopolyyne peak: 
4 h 38 m 38.6 s ,25°35'45", 1950.0) and covering the N=9-8 
rotational transition of C5N. The N=9-8 line is split into 
two fine structure components whose frequencies, derived 
from laboratory measurements, are marked by upward ar- 
rows. The hyperfine structure is too small to be resolved. 
The spectrum was observed by switching the local oscilla- 
tor in frequency by ±200 kHz or ±500 kHz; the position of 
the ghosts of the line in the folded spectrum are indicated 
by downward arrows. 

variations, is estimated to be < 20%. The HC5N and C3N 
line intensities we measure are consistent with those re- 
ported by Tolle et al. (1981) and Guclin et al. (1982). 

Figure 1 (see also Table 1) shows the spectrum covering 
the C5N transitions. We see two 0.3 kms -1 -wide spectral 
lines, each detected at > ha. The lines are separated by 
10.71 ± .01 MHz, which is very close to the value of the 
spin- rotation constant measured by Kasai et al. (1997), 
7 = 10.75 MHz. Their half-power width is comparable to 
the C3N linewidth, 0.23 ± .05 kms -1 . Their rest frequen- 
cies coincide within the small uncertainties with the C5N 
N=9-8 transition frequencies calculated by Kasai et al. 
(1997), if we adopt the source LSR velocity of 5.65 ± 0.05 
kms -1 measured for C3N and C4H (see Guelin et al. 1982). 
Since there are no other comparable lines in the 20 MHz- 
wide spectrum we observed, the probability for a chance 
coincidence is < 10 -6 . We thus conclude that we have 
detected C 5 N in TMC1. 

The 30-m observations were made in Nov. 1997 and 
April 98. The telescope was equipped with two SIS mixer 
receivers with orthogonal polarizations. The zenith atmo- 
spheric opacity was below 0.1 and the system temperature 
T sys = 150 — 200 K. The observations were made by wob- 
bling in azimuth the secondary mirror at a rate of 0.5 Hz 
and with an amplitude of 90". We searched for the N=32- 
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Fig. 2. a): Spectrum observed with the 30-m telescope 
toward IRC+10216. The C 3 H and C 5 N line frequencies 
derived from laboratory measurements are indicated by 
downward arrows. The spectral resolution is 1 MHz (3.3 
kms -1 ). b): The J=4.5-3.5 line emission of C3H in its 
2 n 1/ / 2 ground vibrational state, observed in IRC+10216 
with the IRAM interferometer (Guclin, Lucas & Cer- 
nicharo 1993). The line intensity has been integrated over 
a narrow velocity interval (2 kms -1 ) centred on the star 
velocity. The coordinates represent the offsets in r.a. and 
dec. with respect to the central star. 2c): The central half 
of the spectrum of Fig. 2a, compared to the fitted 3-line 
spectrum (see text). 



31 line, the next two lower lines being partly blended with 
lines of unrelated species. The average spectrum obtained 
toward IRC±10216 is plotted in Figure 2a. It has an r.m.s. 
noise per 1 MHz channel of 0.90 mK in the TjJ scale. 

The carbon-chain molecules observed in IRC±10216 
are concentrated in a thin shell of radius ~ 15" (see 
Fig. 2b). Observed with the 30-m telescope, the profiles 
of their A 3 mm lines have the same cusped shape and 
the same width (29.5 kms -1 ). They are all centred at 
Vlsr= -26.5 kms -1 . The spectrum of Fig. 2a shows two 
cusped lines, near 89.75 GHz and three weaker features 
near 89.80 GHz. The former can be readily assigned to 
the N = 4 -> 3,F = 4.5 -> 3.5 and F = 3.5 -> 2.5 fine 
structure components of C3H in its first excited bending 
state, v 4 = l, 2 £ M (Yamamoto et al. 1990) since they have 
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the right frequencies and the right intensity ratio (1.33±.l, 
whereas the expected LTE ratio is 1.30). 

We have fitted 3 cusped lines (assuming standard line 
shapes and 29.5 kms -1 widths) to the 3 weak spectral 
features of Fig. 2a. The derived intensities and rest fre- 
quencies are given in Table 1; the fitted spectrum is com- 
pared in Fig. 2c to the observed spectrum. The first two 
lines agree in frequency with the (N, J) — (32, 32.5) — > 
(31,31.5) and (32,31.5) -> (31,30.5) transitions of C 5 N 
and very probably arise from this radical. Indeed, although 
the A 3 mm spectrum of IRC+10216 is more crowded than 
the 1 cm spectrum of TMC1, there are not many lines 
that we cannot assign to the rotational transition of a 
known circumstellar molecule: in our 30-m telescope spec- 
tral survey of IRC+10216, the density of unidentifed lines 
stronger than > 3 mK is only of 1 per 80 MHz (see Guelin 
et al. 1997). The probability that two unrelated lines lie 
within 1 MHz from the C5N lines is thus lower than 1 per 
one thousand. Whereas the first two features can be ten- 
tatively assigned to C5N, the third, which is weaker and 
lies 7 MHz higher in frequency, remains unidentified. 

3. The abundance of C 5 N 

In the direction of the TMC 1 cyanopolyyne peak, the ro- 
tational populations of HC5N and HC3N can be described 
by Boltzmann distributions with rotation temperatures of 
7 - 10 K (see e.g. Takano et al. 1997). We adopt there- 
fore Trot — 8 K. The H 2 column density in this direction is 
N(H 2 )= 10 22 cm" 2 (Cernicharo & Guelin 1987). That of 
C5N can be calculated from the standard expression for 
the optically thin lines of thermalizcd linear molecules: 

N(C 5 N) = 0.710 16 -p^% M-W^S / T MB dv, 

1 rot l-bg J 

where 

T'=(hv/k)(e& -l)- 1 , 

and where E fT M Bdv = 0.014 Kkms -1 is the sum of 
the integrated intensities of the two doublet components. 
In equation [1], N is in cm -2 , the dipole moment /xq = 
3.385 D (Botschwina 1996) in debye, and the line fre- 
quency v — 25.25 GHz in gigahertz. We find: 
N(C 5 N) = 3.110 n cm- 2 , 
£(C 5 N) = N(C 5 N)/N(H 2 ) ~ 310- 11 . 

The abundance of C5N can be compared to those of the 
related species HC5N and C3N in the light of the chemical 
model predictions. We calculate first the abundance of the 
13 C isotopomers of HC5N, whose 1 cm lines are optically 
thin, by setting in equation [1] fi = 4.33 D, v = 23.7 GHz 
and J TmbC^v = 0.065 Kkms -1 , which is the average of the 
integrated intensities of the J=9-8 lines of HC 13 CC3N and 
HC4 13 CN, observed with the Effelsberg telescope (Takano 
et al. 1998). We find N(HC 4 13 CN)= 1.0 10 12 cur 2 , from 
which we derive N(HC 5 N)= 7 10 13 cm" 2 , adopting the 



'standard' elemental abundance ratio 12 C/ 13 C~70 in the 
local interstellar medium (see Wilson & Rood 1994). The 
value of 7 10 13 cm~ 2 is close to the value of Suzuki et al. 
(1992), as well as to the value we estimate with an LVG 
code from the main isotopomer J = 9 — > 8 line intensity 
(5 10 13 cm~ 2 ). We arrive at: 
N(HC 4 13 CN/C 5 N)= 3 
N(HC 5 N/C 5 N)~ 200. 

For C 3 N, we find, using in eq. [1] T rot = 8K, /j, = 2.84 
D (Pauzat et al. 1991), and the line parameters of Table 
1: N(C 3 N)= 8.2 10 12 cm" 2 . Finally, we take for HC 3 N the 
column densities derived by Takano et al. (1998), 
N(HC 13 CCN)= 2.1 10 12 cm" 2 and N(H 12 C 3 N) = 1.6 10 14 
cm -2 . This yields: 
N(HC 13 CCN/C 3 N)= 0.26 
N(H 12 C 3 N/C 3 N)= 19. 

The C5N/HC5N abundance ratio is an order of magnitude 
smaller than the C 3 N/HC 3 N ratio. According to neutral- 
neutral gas phase chemical models, the cyanopolyyne 
mole-cules are mainly formed by reactions of N and CN 
with polyacetylenes or polyacetylic ions (Herbst & Leung 
1990). The cyanopolyyne radicals are formed in TMC1 by 
the reaction of atomic C with cyanopolyynes (Herbst et 
al. 1994), and in IRC+10216 by photodissociation (Cher- 
chneff & Glassgold 1993). They are destroyed by reactions 
with N atoms, O atoms, polyacetylenes and photodissoci- 
ation. 

According to model predictions, HC3N/HC5N varies in 
TMC1 by several orders of magnitude between the "early 
times" and steady state. The C3N/HC3N and C 5 N/HC 5 N 
ratios, on the other hand, remain constant within a factor 
of 2. They are comprised between 0.1 and 0.2 (see Table 3 
of Herbst et al. 1994). Whereas the observed C 3 N/HC 3 N 
ratio (0.2, Cernicharo et al. 1987) agrees with the pre- 
dicted one, C5N/HC5N is more than one order of mag- 
nitude too low. The most recent models, which take into 
account the destruction of acetylene and polyacetylenes 
by C atoms and the reaction of CN with O, form enough 
C 5 N but too little C 3 N, HC 3 N and HC 5 N (Herbst et al. 
1994). 

In the case of IRC+10216, we adopt for C 5 N the 
same rotation temperature as for HC5N (T rot = 29 K, 
Kawaguchi et al. 1996). We then derive from the inte- 
grated intensities of Table 1 a line-of-sight column den- 
sity in the direction of the central star (twice the ra- 
dial column density across the shell) N(CsN)= 6 10 12 
cm -2 . This is ~ 50 times less than the column densi- 
ties of C 3 N and HC5N. Here also, C5N is underabundant 
with respect to model predictions (Cherchneff & Glassgold 
1993), and C5N/HC5N one order of magnitude smaller 
than C 3 N/HC 3 N. 

The unexpectedly low C5N/HC5N abundance ratio 
in found both sources shows that the formation of long 
carbon-chain molecules is not fully understood, and that is 
difficult to predict the abundances of unobserved species. 
The very long carbon chains could be more abundant than 
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Table 1: Observed line parameters 



Rest, frequ. 


Obs.-Calc. 


Transition 


Species 


J T M sdv 


(MHz) 


(MHz) 


7V,J-> N',J' 


(mK.kms- 1 ) 


TMC1 


19799.956 


.005 


2,1.5,1-5-1,0.5,0.5 


C 3 N 


23 (8) 


19800.121 


.000 


2,1.5,2.5-1.0.5,1.5 


C 3 N 


58 (8) 


19780.801 


.001 


2,2.5,2.5-1,1.5,1.5 


C 3 N 


60 (8) 


19780.826 


.000 


2,2.5,1.5-1,1.5,0.5 


C 3 N 


34.5 (8) 


19781.096 


.002 


2,2.5,5.5-1,1.5,1.5 


C 3 N 


86 (8) 


23963.897 


0.000 


9-8 


HC 5 N 


2240 (100) 


25249.938 (4) 


.018 (40) 


9,9.5-8,8.5 


C 5 N 


7.3 (9) 


25260.649 (4) 


-.017 (40) 


9,8.5-8,7.5 


C 5 N 


6.4 (9) 


IRC+10216 


89730.54(10) 


-.07 (10)" 


4,4.5 - 3,3.5 


C 3 H(i^E") 


880(20) 


89759.17(12) 


-.18 (12) 6 


4,3.5 - 3,2.5 


CsH^E") 


660(20) 


89785.6 (4) 


.4 (13) 


32,32.5-31,31.5 


C 5 N 


95 (15) 


89797.0 (3) 


1.1(13) 


32,31.5-31,30.5 


C 5 N 


105(20) 


89804.0 (10) 






U 


75 (20) 



Notes to the Table: 

b : weighted average of the two blended hyperfine components. The calculated frequencies are taken from Kasai et al. (1997) 
for C 5 N, Yamamoto et al. (1990) for C 3 H, and Guelin et al. (1982) for C 3 N (sec also Gottlieb et al. 1983). The observed 
rest frequencies of C 3 N and C5N were derived assuming Vlsr= 5.65 kms -1 in TMC 1 and -26.5 kms" 1 in IRC+10216; 
that of HC5N was derived assuming Vlsr= 5.75 kms -1 (see Guelin et al. 1982). The number in parenthesis represent 
the r.m.s. uncertainty on the last digit. The uncertainties on the line integrated intensities, given in the Table, do not 
include the calibration uncertainty which is 10 — 20%. Note that the intensity ratio between the rotational transitions in the 
C 3 H ground state and in the 2 E M excited bending state, hence presumably the population ratio between these two states, is ~ 10. 



expected. We note, however, that the species of medium 
size, such as the chains consisting of 4 — 7 C,N, or O atoms 
and the rings with less than 10 heavy atoms, which would 
give rise in TMC 1 to a rich centimetric spectrum, are 
probably not very abundant in that source: we have cov- 
ered so far a 100 MHz-wide band in TMC 1 with a very 
good sensitivity and did not detect any unidentified line 
down to a level of 10 mK. The acetylenic chains and cu- 
mulene carbenes appear in this respect exceptional. 
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